"Modified" BCS theory at finite temperature and its flaws 
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Drawbacks of a recently suggested modified BCS theory at finite temperature in [Phys. Rev. C 
67, Of 4304 (2003)] are discussed. We conclude that the thermal behavior of the pairing correlations 
is not under control within this model. 

PACS numbers: 21.60.-n, 24.10.Pa 



A modified BCS (MBCS) model has been suggested recently pj to improve the BCS theory at finite temperature 
for nuclei. It has been claimed that the MBCS "washes out the sharp superfluid-normal phase transition" and u the 
fluctuations of particle number also become more suppressed especially at high temperature" . In this Comment we 
argue that the corresponding equations cannot be treated as an extension of the modified BCS approach to finite 
temperature. Consequently, the "physical effects" reported in this article are only the results of incorrect physical 
treatment and inconsistencies. 

The MBCS equations at T ^ are obtained applying a temperature-dependent unitary transformation to the 
Bogoliubov quasiparticles {a] m , ctj m }, thus transforming them into the new bar quasiparticles {a] m , a.j m }: 

*) m = V 71 - njajm + (-l) J '~ m ^<*i-m, (1) 

where rij arc the thermal Fermi-Dirac occupation numbers for the Bogoliubov quasiparticles. The authors claim: 
"The MBCS equations . . . has exactly the same form as the standard BCS equations (At=o = G^2j ^j u j v j) where 
the coefficients Uj and Vj are replaced with Uj and Vj ": 



A(T) = (: Y, U - 



(1 - 2rij)ujVj 

(2) 



Although formally the conventional BCS (T = 0) and MBCS equations look similar, there is one essential difference 
between them. As far as Uj and Vj coefficients are positively defined, the BCS gap A = (breaking of the Cooper 
pairing) corresponds to the trivial solution {uj, Vj} — {0(1), 1(0)} for all j when all levels give zero contribution to 
the gap. It is not difficult to show analytically that the MBCS gap equation J2J does not have such a trivial solution 
{Uj, Vj} = {0(1), 1(0)}. Indeed, from Eq. (38) in Ref. 



uj = Uj v/1 - rij + vj y/nj, Vj = Vj s/l-rij - Uj \fnj (3) 

the trivial solution {u,j, Vj} — {0(1), 1(0)} corresponds to: 

Uj = */ 1 — nj ; Vj = \frtj - particles 

y/1 — nj - holes 



u 3 = -V n i> v 



and contradicts the positive definition of Uj . 

One may notice from Eq. that Vj coefficients become negative above a certain temperature for particle levels 
since Vj << Uj for them. Then the MBCS pairing gap A receives positive contribution from hole levels and negative 
contribution from particle levels. This leads to quite a strange thermal behavior of A(T). 

On page 5 of Ref. Q we find a discussion of the MBCS critical temperature T c defined as A(T C ) = 0. According 
to the authors this temperature is reached when the "new" (second) term in Eq. @ is equal to the conventional 
one. Numeric calculations (see below) show that two terms in Eq. J2J) compensate each other around the critical 
temperature of the conventional BCS T c w 0.57 • At=o for particle levels and they never do for hole levels because 
the second term in Eq. J2J is always negative for holes. 

As clear from the above paragraphs, A vanishes at T c only when a negative contribution from particles and positive 
contribution from holes cancel each other (notice the difference with the conventional BCS A = 0). However if this 
happens, at higher T the balance appears to be broken and A becomes finite again. 

To be more specific, we plot in Fig.^, the MBCS pairing gap as function of temperature in 76 Ni One notices 
that it reaches zero at T c w 2.1 MeV and continues to decrease. The authors of Ref. define an unphysical negative 
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FIG. 1: (a) the MBCS pairing gap; (b) Vj and Vj coefficients for particle levels near the Fermi surface: 2d 5 / 2 - solid lines, 3si/2 
- dashed lines; (c) two functions for particle-number fluctuations v/ SNf in Eq. EJ • Critical temperates T c and T c are indicated 
by vertical lines. See text for details. 



gap of the MBCS as the results "no longer reliable" . Since the authors never discuss any physical limitations of their 
model, it is not clear upon what basis this conclusion is made. In Fig.^D we present the Vj and Vj coefficients for two 
particle (2d 5 / 2 j 3s!/ 2 ) levels near the Fermi surface. They become negative around T c . However, one may notice that 
nothing special occurs at f c , thus confirming the "unreliability" of the method above the "new" critical temperature. 

The problem with a negative gap can be removed by increasing the role of hole levels. It may be achieved, e.g., 
by including the levels from the shell N=0-28 in a single-particle spectrum in addition to the levels near the Fermi 
surface accounted for in Ref. 

(with a proper renormalization of the pairing strength). Then, another MBCS crack 
appears: the gap starts to continously increase above a certain temperature remaining always positive. Similar effect 
in 84 Ni is shown in Figs. 2h and 3 of Ref. But we find it impossible for the MBCS gap to reach zero and remain 
zero at higher temperature, thus, indicating the superfluid-normal phase transition. 

In view of all the above, the author's conclusions about superfluid-normal phase transition within the MBCS makes 
no sense from our point of view. 

The other main goal of the MBCS is to reduce the particle number fluctuations in the ground state. The authors 
show in Fig. 5 (solid line) of Ref. Q that this quantity 

^ 2 = 4^^=A 2 ^^. (4) 

3 3 

is reducing with the temperature and apparently vanishes. The 5N 2 has been calculated from the r.h.s. part of 
Eq. © (or Eq.(41) in Ref. Q). It remains unclear how transformation in Eq. Q has been obtained within the MBCS 
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at T ^ 0. We have calculated both the left and right functions. Numeric results presented in Fig. ^ reveal their 
rather different behavior. It means that the conclusion on suppression of the particle number fluctuations with a 
temperature within the MBCS in Ref. pj is based on calculations from the equation which is not valid. 

A legitimate question is what is wrong with the MBCS. In Sec. II. D and Appendix A of Ref. [l| the well-known 
statistical approach in a many-body theory at finite temperature is sketched. However, deriving the MBCS equations 
the authors of Ref. [jj do not use this method. Instead, they formally follow the scheme formulated by them for zero 
temperature. Namely, they find a minimum of the expectation value of the BCS Hamiltonian over the temperature 
dependent vacuum state |0). They imply for reasons unclear that selecting rij(T) in the form of the quasiparticle 
thermal occupation numbers is equivalent to averaging over the grand canonical ensemble (as one reads on page 5 of 
Ref. Q). 

We find it difficult to understand: the authors are aware that a quantum mechanical ground state, for which 
expectation value of the Hamiltonian equals the statistical average over the grand canonical ensemble, "does not exist 
in the physical space spanned by eigen vectors" of the pairing Hamiltonian and why it is so (see footnote on page 3). 
Nevertheless, the basic starting point of the MBCS at T ^ is certainly to introduce such a vacuum state 0) in 
Eq. (33). 

This means that the MBCS cannot be treated as a theory with the thermal behavior of the pairing correlations 
under control. It is not surprising that the MBCS equations (Eqs. (39-40) and (43-44) in Ref.|l() differ from the 
conventional ones by extra terms but their origin is only inconsistency in treatment of the physical object. 

To conclude, we doubt that the modified BCS in Ref. [l| with the above discussed flaws is a step forward in the 
developing of the conventional BCS at finite temperature. 
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[2] Our code reproduces excellently all the results in Ref. Q for all isotopes. We present them for 76 Ni only (from the resonant 
continuum MBCS calculations). The ones in 68 " 76 Ni and 80 - 82 Ni are qualitatively very similar and lead to the same 
conclusions. 



